Two qubits can be entangled in two distinct temperature regions 
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We have found that for a wide range of two-qubit Hamiltonians the canonical-ensemble thermal 
state is entangled in two distinct temperature regions. In most cases the ground state is entangled; 
however we have also found an example where the ground state is separable and there are still two 
regions. This demonstrates that the qualitative behavior of entanglement with temperature can 
be much more complicated than might otherwise have been expected; it is not simply determined 
by the entanglement of the ground state, even for the simple case of two qubits. Furthermore, we 
prove a finite bound on the number of possible entangled regions for two qubits, thus showing that 
arbitrarily many transitions from entanglement to separability are not possible. We also provide 
an elementary proof that the spectrum of the thermal state at a lower temperature majorizes that 
at a higher temperature, for any Hamiltonian, and use this result to show that only one entangled 
region is possible for the special case of Hamiltonians without magnetic fields. 

PACS numbers: 03.67.Mn, 03.65.Ud 



I. INTRODUCTION 

A topic that has emerged recently within the field of 
quantum information is the study of entanglement in spin 
systems (see [J, |2j for early examples). Entanglement 
is often considered necessarily a low-temperature phe- 
nomenon that becomes less important as the temperature 
is increased. It was therefore surprising that in [lj an ex- 
ample for two qubits was given where the ground state 
(T = 0) is separable but the thermal state is entangled at 
higher temperatures. This behavior can be understood as 
due to the presence of low-lying excited states that are en- 
tangled. Thus at least two qualitatively-different entan- 
glement scenarios are possible for two qubits (apart from 
the uninteresting case of no entanglement at any temper- 
ature); (i) the ground state is entangled, and hence the 
thermal state is entangled at low temperatures up to a 
critical temperature, T5, above which it is separable, and 
(ii) the ground state is separable, but the thermal state 
is entangled for temperatures within some finite range 
(Ti, Ts) (and separable again above Tg). 

The question we address here is what other entan- 
glement scenarios are possible for two qubits. We were 
stimulated to ask this question by which studies the 
generic behavior of thermal entanglement as a function of 
temperature. There it was shown that the generic behav- 
ior is closed intervals in temperature where the thermal 
state is separable interspersed with open intervals of en- 
tanglement. Examples of the two types mentioned above 
were given for two qubits. Also two examples of qubit- 
qutrit systems {d\ — 2 and d 2 — 3, where dj, j = 1,2 
are the dimensions of the two subsystems) were given 
where there are two distinct entangled regions; in one 
case the ground state is entangled, and in the other case 
the ground state is separable. 
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Non-monotonic behavior of thermal entanglement has 
also been observed for qubit spin chains Ref. Q 

studies the reduced state of two nearby qubits in particu- 
lar spin chains and finds examples of two distinct entan- 
gled regions in temperature. Ref. Q uses a multipartite 
entanglement measure, and observes three regions. In Q 
transitions from separability to entanglement are studied 
as a type of phase transition using geometric arguments 
about the set of separable states. 

Here we restrict to the case of just two qubits (as op- 
posed to the reduced state of two qubits out of a large 
chain), and present a class of Hamiltonians for which 
most have a value of the magnetic field strength such 
that the ground state is entangled and there are two 
entangled regions. In addition, we present an example 
of a Hamiltonian outside this class that has a separable 
ground state and, again, two entangled regions. Thus we 
find that all the classes of behavior for the thermal states 
of qubit-qutrit systems found in Ref. Q are also observed 
for the two-qubit case. 

These results raise the question of how many distinct 
"entangled regions" in temperature are possible. Our 
numerical search failed to find Hamiltonians with more 
entangled regions for two qubits, indicating that this is 
the most complicated behavior. We show that for a class 
of commonly considered Hamiltonians — those without 
magnetic fields — it is impossible to obtain more than 
one entangled region. In addition we derive upper bounds 
on the number of entangled regions in the general case. 

This paper is set out as follows. In Sec. [H] we present 
results for the dimer case of the spin-chain Hamiltonian 
studied in Ref. Q, then give our general class of two- 
qubit Hamiltonians. In Sec. IIIII we give an example of a 
Hamiltonian in our class that does not exhibit two en- 
tangled regions, and show that small perturbations are 
sufficient to give two entangled regions. In Sec. IIVI we 
give our example with a separable ground state and two 
entangled regions. We derive bounds on the total num- 
ber of entangled regions possible for two qubits in Sec.lVl 
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specialize to the case with zero magnetic field in Sec. I VII 
and then conclude in Sec. IV11I 

II. GENERAL CASE 

The two-qubit case of the XYX spin-chain studied in 
Ref. is of the following form: 

H = -J[X 1 X 2 - Z 1 Z 2 + h{Z 1 + Z 2 )]. (1) 

where {Xj,Yj, Zj} are the Pauli sigma matrices acting 
on qubits j — 1,2, J is a coupling constant with di- 
mensions of energy and h is a dimensionless parameter 
corresponding to the magnetic fields experienced by the 
qubits. The results for this Hamiltonian on two-qubits 
were given in the inset of Fig. 4 in an early version of this 
paper • This figure showed that two entangled regions 
were obtained, though this aspect of the results was not 
discussed. 

An isolated system (i.e. not exchanging particles with 
the environment) in thermal equilibrium with a bath at 
temperature T will reach the canonical-ensemble thermal 
state p(T) given by 

P(T) = (2) 

where j3 = 1/fegT and Z = Tr[exp(— 0H)] is the parti- 
tion function. 

In Fig. n we plot the concurrence § as a function 
of (scaled) temperature and the scaled magnetic field h 
for the canonical-ensemble thermal state Hi . We see 
that there is a significant range of values of h (approxi- 
mately \J~2 < h < 2.36) such that the entanglement as a 
function of temperature behaves as claimed — there are 
two entangled regions, one at low temperatures, [0,Ti), 
and then another distinct region at higher temperatures, 
(T 2t T s ). 

A general two-qubit Hamiltonian may be written in 
the form (omitting any global phase shift) 

H = X"[RX 2 + AX X + BX 2 , (3) 

where R is a real 3x3 matrix, A and B are real row 
vectors, and Xf, = [Xf, Yk Zk] T , k G {1 7 2}. It is possi- 
ble to add local unitary operations before and after the 
Hamiltonian without affecting the entanglement of the 
thermal state. This is because 

e -pvi*v a Hv}9v2 =Vi9 v 2e -PH V i ® vl ( 4 ) 

where V\ and V 2 are local unitary operations on subsys- 
tems 1 and 2. These local unitaries act to rotate the 
Xk vectors, giving VkXkV^ = OkXk for some orthogonal 
matrices Ok- Thus 

H' = V x ® V 2 HVl <g> V$ 

= XjOlR0 2 X 2 + AOxXx + B0 2 X 2 . (5) 
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FIG. 1: Concurrence as a function of temperature and h for 
two qubits coupled according to 

If the local operations are chosen such that 0\ and 2 
are the orthogonal matrices which result from a singular- 
value decomposition of i?, then 0\R0 2 is a positive di- 
agonal matrix 

Hence W is of the form 

H' = J[a x X x X 2 + a y Y x Y 2 + a z Z,Z 2 + h(f3 xl X 1 
+ P x2 X 2 + f3 yX Y x + l3 y2 Y 2 + zl Z, + I3 Z2 Z 2 % (6) 

where the aj are positive, and the ftjk are real. Note that 
the local unitaries do not remove the local component 
of the Hamiltonian. If it were possible to use different 
local unitaries before and after the Hamiltonian, the local 
component of the Hamiltonian could be removed entirely 
[Tfl |. However, this would change the entanglement of the 
thermal state. 

Hamiltonians of the form (JSJ are the most general two- 
qubit Hamiltonians for the problem of thermal entangle- 
ment. Now we introduce a class of Hamiltonians which 
is slightly restricted, in that we require (3j\ = (3j 2 = f3j. 
This is equivalent to requiring the magnetic field to be 
homogeneous. These Hamiltonians may be written as 

H = Jia^X-2 + a v Y x Yi + a z Z x Z 2 

+h\fl x {X x + X 2 ) + 0y(Yi + Y 2 ) + (5 Z {Z X + Z 2 )}} . 

(7) 

As before, the aj are positive, and the 0j are real. To de- 
termine properties of these Hamiltonians, random Hamil- 
tonians were generated, and for each it was determined if 
there exists a value of h such that there are two entangled 
regions. The aj were chosen at random in the interval 
[0, 1), and the (3j at random in the interval (—1, 1). From 
a sample of 2 20 of these Hamiltonians, it was found that 
all had a value of h such that there are two entangled 
regions. 

Arbitrary two-qubit Hamiltonians were also tested. 
These were generated according to the Gaussian unitary 



3 



20 r 

18 
16 
14 
12 
10 

8 

6 

4- 

2 



0.5 



k B T/J 



1.5 



2.5 



10.8 



10.7 




FIG. 2: Concurrence as a function of temperature and h for 
two qubits coupled according to @ with 7 = 10 -6 . 



FIG. 3: Concurrence as a function of temperature and h for 
two qubits coupled according to (I1UI) with 8 = 10~ 6 . 



ensemble. Each Hamiltonian was divided into a local 
part Hl and a nonlocal part if/v, and Hamiltonians of 
the form = + hH^ were tested. It was found 
that, out of 2 12 samples, there were 106 such that there 
was a value of h for which Hh has two entangled regions. 
This gives the overall probability for this behavior for 
arbitrary Hamiltonians as 2.59 ± 0.25%. 



For 8 <C 1 this is the XY interaction with a misaligned 
transverse magnetic field. The concurrence for 8 = 10~ 6 
is shown in Fig. |3| Even with this very small misalign- 
ment in the magnetic field, the distinct entangled regions 
are again seen. 



III. EXAMPLES 

Although we have shown that it is extremely common 
for Hamiltonians of the form J7J to have two entangled 
regions, not all exhibit this behavior. For example, con- 
sider the XY interaction with a magnetic field, as studied 
by Wang \u&: 

H = J{X 1 X 2 + Y 1 Y 2 + h(Z 1 + Z 2 )}}. (8) 

Wang found that it was possible for the thermal entangle- 
ment to be zero for T = but nonzero for T > 0. Wang 
also considered the anisotropic XY interaction, but with- 
out a magnetic field. In neither case were two regions 
found. 

It turns out that we can vary the Hamiltonians very 
slightly from this example, and again recover the two 
entangled regions. For example, consider the anisotropic 
XY interaction 

H = J [(1 + 7)^1^2 + (1 - l)YxY 2 + h{Z x + Z 2 )\ . (9) 

For 7 = this is the Hamiltonian of Eq. © ■ However, for 
7 equal to just 10 -6 , we again recover the two entangled 
regions (see Fig- EJ- 

Another perturbation which recovers the two entan- 
gled regions is that where the magnetic field is not ex- 
actly aligned on the z-axis: 

H = J {X X X 2 + Y X Y 2 + h[Z x +Z 2 + 8{X 1 + X 2 )}} . 

(10) 



IV. SEPARABLE GROUND STATE 

The next most complicated case is that where the 
ground state is separable, so the thermal state is separa- 
ble at T = 0, but there are still two entangled regions. As 
local unitaries do not alter the entanglement, one can ar- 
bitrarily choose the separable ground state without loss 
of generality. Thus to numerically search for such ex- 
amples we took the ground state to be |00). The other 
eigenstates and eigenenergies were then chosen at ran- 
dom. 

The example found was (after rounding the coeffi- 
cients) 

H = 0.006(A!A 2 + YiY 2 ) + 0.03(A!y 2 - YiX 2 ) 

+ nm{z x x 2 - x 2 ) + {z x y 2 - r 2 )/io 

+ (X 1 Z 2 - Xi)/U + Z ± Z 2 /7 - Zi/4 - Z a /5. (11) 

The concurrence as a function of temperature is as shown 
in Fig.0] There are two distinct regions of entanglement, 
with a separable ground state. Note that there appears 
to be a finite region without entanglement for low tem- 
perature. However, for much of this region the concur- 
rence is extremely small (less than 10 -20 ) but nonzero. 
This indicates that the thermal state may be completely 
separable only for zero temperature, and the entangle- 
ment for small temperatures is not observed due to finite 
precision. 
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FIG. 4: Concurrence as a function of temperature for two- 
qubit thermal state under the Hamiltonian li lt . 



V. BOUND ON NUMBER OF ENTANGLED 
REGIONS 

We used numerical techniques to search for examples of 
more complicated scenarios from the hierarchy, i.e. three 
or more entangled regions, but were unable to find any. 
Of course no numerical technique can be exhaustive so it 
remains an intriguing possibility that even more entan- 
gled regions are possible, even for two qubits. We now 
show that there is, in fact, a finite upper bound, 17, on 
the number of entangled regions for two qubits. How- 
ever it remains entirely plausible that this bound is not 
tight and the above examples of two entangled regions 
represent the most complicated behavior possible for two 
qubits. 

To analytically bound the number of entangled regions 
for two qubits we use the well-known fact that a two- 
qubit mixed state is entangled or separable depending 
on whether its partial transpose with respect to one of 
the qubits has a negative eigenvalue or not. Therefore, 
by solving 

det[p{T) TA ] = 0, (12) 

where Ta denotes the partial transpose with respect to 
the first subsystem, we find the transitions between en- 
tangled and separable thermal states. We may scale the 
energies so that the minimum energy eigenvalue is zero. 
If we multiply by Z, the resulting equation is polynomial 
in er x l kaT with noninteger powers and 35 terms [TtI ]. 

Provided the ratios of the energy eigenvalues are all ra- 
tional, the polynomial has integer powers in x = e~ r / kBT 
for some constant r. To place a bound on the number of 
solutions, we first take the derivative of the polynomial, 
then apply Descartes' rule of signs. The derivative of the 
polynomial has no more than 35 terms, and so has no 
more than 34 sign changes. By Descartes' rule there are 



no more than 34 (positive) zeros of the derivative, and 
no more than 34 turning points of the polynomial. 

It is easily seen that, provided the derivative has no 
more than 34 zeros, there are no more than 17 regions 
where the polynomial is negative. In the two-qubit case, 
the partial transpose has no more than one negative 
eigenvalue thus the determinant is negative if and 
only if the partial transpose is negative. Hence there can 
be no more than 17 entangled regions. Thus we obtain a 
finite limit on the number of entangled regions, though 
this is much larger than the number of intervals which 
have been found numerically. In practice the number of 
sign changes is likely to be far less than 34, although we 
do not see a way of showing this analytically. 

In the case where there are irrational ratios of the en- 
ergy eigenvalues, the situation is more complicated be- 
cause the polynomial has noninteger powers. In this case, 
we can achieve an arbitrarily close approximation of the 
Hamiltonian with rational energies. In the case where a 
function / is the limit of a sequence of functions /„, it 
is not possible for / to have more turning points than 
/„. The only situation where /' (where the prime indi- 
cates the derivative) can have more zeros than f n is when 
f' n has an extremum which approaches zero in the limit 
n — > co, and is only exactly equal to zero for /'. However, 
this zero would correspond to a point of inflection, rather 
than an extremum, for /. Thus we find that the polyno- 
mial in e~ 1 / kBT can have no more than 17 regions where 
it is negative, and the limit on the number of entangled 
regions must hold for irrational powers also. 

For a qubit coupled to a qutrit entangled mixed states 
must still have a non-positive partial transpose |l3j | , how- 
ever there is a complication due to the fact that the par- 
tial transpose can have more than one negative eigen- 
value. The main problem in this case is that, at a point 
where det[p(T) TA ] changes sign, the state could be sep- 
arable for det[p(T) TA ] — 0, but entangled for slightly 
higher or lower temperatures. This could happen if one 
of the eigenvalues passes from positive to negative, while 
another passes from negative to zero to negative. 

However, despite this possibility it can be seen that 
the number of turning points of det[p(T) TA ] still pro- 
vides an upper bound on the number of entangled re- 
gions. For a qubit-qutrit system there are six energy 
levels and so up to 462 terms in the polynomial corre- 
sponding to det[p(T) TA ] (assuming rational eigenvalues). 
There are therefore no more than 461 turning points, 
even in the limit of irrational eigenvalues. Combined 
with the fact that the state must be separable at high 
temperature, this implies that the number of entangled 
regions can be no higher than 462 (some of them may 
be separated by single points in temperature where the 
system is separable). 

More generally, for the case of two subsystems of ar- 
bitrary dimensions, d\ and c?2, one might hope to put a 
finite upper bound on the number of entangled regions as 
a function of d\ and d^. However, in higher dimensions 
entangled mixed states do not necessarily have a non- 
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positive partial transpose. In fact it has recently been 
shown that even the problem of distinguishing separable 
and entangled mixed states is TVP-hard in arbitrary di- 
mension |14|. It is therefore unlikely that this approach 
will yield upper bounds for higher dimensional systems. 



VI. ZERO MAGNETIC FIELD 

Although we were unable to definitively answer the 
question of what entanglement scenarios are possible for 
an arbitrary two-qubit Hamiltonian, we can for a certain 
class of Hamiltonian — those that have no local terms 
(corresponding to a magnetic field), and only interaction 
terms. A Hamiltonian without local terms may be writ- 
ten in the form 



Hence 



H — X 1 RX 2 , 



(13) 



As in Sec.[nl we can apply local unitaries without altering 
the entanglement of the thermal state. These simplify the 
Hamiltonian to a form that is diagonal in the Bell Basis. 

H' = V X ®V 2 HV}®VZ = J{a x X 1 X 2 +a v Y 1 Y 2 +a z Z 1 Z 2 ). 

(14) 

The Bell basis is a set of maximally entangled states 



1^) = ^ CiT>]T>±U>|l» , 
1^) = -^(|T)U>±U>|T» , 



(15) 
(16) 



where Z\\) — ||), Z\{) — —\\) (up and down spins if the 
qubit is a spin 1/2 quantum system). 

To determine the behavior of the entanglement, we 
compare the state at two different temperatures, Ti and 
T 2 , such that T\ < T 2 , We first note that the eigen- 
values for p(T 2 ) are majorized by those for p{T\). It is 
straightforward to show this result for all bipartite ther- 
mal states (not just those for two-qubit systems). For 
Ti < T 2 , /3i > (5 2 . Therefore, for AE > (<) 0, we have 
e -/3 2 AB > e -p 1 AE _ Taking the energy eigenvalues Ei 
to be sorted into nondescending order, we have 



E^ 1 



(Ei-E k ) 



> 



E< 



-/3 2 (Ei-E k 



(17) 



e -fr(Ei-E h ) > e -^- E "\ (18) 



i=Jt+l 
Multiplying gives 

k d 



i=fc+l 



, (19) 



=fe+i 
d 



i=l 
k 



i=k+l 



Y e- piE ' Yl e ~ P2Ei ^ E e ~ P2El E e ^ lEi ■ ( 2 °) 



(21) 



which is the result claimed. 

Now, for density operators p\ and p 2 such that the 
eigenvalues for p 2 are majorized by those for pi, we have 



P2 = Y,Pi U }pi U h 



(22) 



where the unitaries Uj permute the eigenstates |18| . 

For the specific case where the Hamiltonian is diagonal 
in the Bell basis, the eigenstates are just the Bell basis 
states. In this case, if the state p\ is separable, then 
all states Ujp\Uj obtained by permuting the eigenstates 
are also separable. To show this, we first note that it 
is not necessary to preserve phase when permuting the 
eigenstates, because any phase cancels out in the density 
matrix. In order to permute the eigenstates (without re- 
gard for phase) , it is sufficient to show that it is possible 
to perform three swaps between eigenstates. All permu- 
tations may be constructed from these three swaps. We 
may obtain three swaps between Bell basis pairs using 
local unitaries as follows: 



r 



e i7rZ/4 <g) e M 
e-" z / 4 (g) e 



Z/4 



(23) 
(24) 
(25) 



where H — (X+Z)/y/2 is the Hadamard operation. Thus 
we find that it is possible to perform any permutation of 
the Bell basis using local operations, so if p\ is separable, 
each of the states l/JpiUj is separable. Hence the state 
p 2 must be separable. 

Thus, for Hamiltonians that are diagonal in the Bell 
basis, the eigenvalues of p(T 2 ) for T 2 > T\ are majorized 
by those for p(T\), so if p(T\) is separable, then so is 
p(T 2 ). Thus we can not have a situation where the ther- 
mal entanglement increases with temperature. As we 
may simplify any two-qubit Hamiltonian without local 
terms to a form which is diagonal in the Bell basis, this 
result holds for all two-qubit Hamiltonians without local 
terms. 



VII. CONCLUSIONS 

We have presented a class of Hamiltonians for which 
almost all examples have a value of the magnetic field 
such that the thermal state has two distinct entangled 
regions. One example of this class previously appeared 
as a figure in an online paper 7] , but this aspect of the 
results was not discussed explicitly. This result is some- 
what surprising as one may have expected that the small 
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Hilbert space for two qubits would mean that only one 
entangled region were possible. 

There are, however, particular cases from this class 
where distinct regions do not occur, for example the 
isotropic XY interaction with a transverse magnetic field 
However, we find that if the interaction is perturbed 
only slightly by making it anisotropic or misaligning the 
magnetic field, distinct regions do occur. This suggests 
that those cases where the distinct regions do not occur 
are a set of measure zero in this class. 

It is also possible for there to be two entangled regions 
when the ground state is separable. In contrast to the 
case where there are two regions and an entangled ground 
state, this behavior is extremely rare. It was necessary to 
test millions of Hamiltonians before an example of this 
form was found. 

We have also shown that certain features of the ex- 
amples are necessary in order to observe the distinct en- 
tangled regions. We proved that for Hamiltonians with- 
out local terms (i.e. no magnetic field) the entanglement 



must necessarily decrease with increasing temperature, 
so only one entangled region is possible (at low temper- 
atures). 

For general two-qubit Hamiltonians we showed, by con- 
sidering zeros of the determinant of the partial transpose 
of the thermal density matrix, that there can be no more 
than 17 entangled regions. Thus arbitrarily many transi- 
tions from entanglement to separability are not possible 
for two qubits (or a qubit and a qutrit). 
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